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ABSTRACT 

We compute grids of radiative-convective model atmospheres for Jupiter, Saturn, 
Uranus, and Neptune over a range of intrinsic fluxes and surface gravities. The atmo- 
sphere grids serve as an upper boundary condition for models of the thermal evolution 
of the planets. Unlike previous work, we customize these grids for the specific proper- 
ties of each planet, including the appropriate chemical abundances and incident fluxes 
as a function of solar system age. Using these grids, we compute new models of the 
thermal evolution of the major planets in an attempt to match their measured luminosi- 
ties at their known ages. Compared to previous work, we find longer cooling times, 
predominantly due to higher atmospheric opacity at young ages. For all planets, we 
employ simple "standard" cooling models that feature adiabatic temperature gradients 
in the interior H/He and water layers, and an initially hot starting point for the calcula- 
tion of subsequent cooling. For Jupiter we find a model cooling age ~10% longer than 
previous work, a modest quantitative difference. This may indicate that the hydrogen 
equation of state used here overestimates the temperatures in the deep interior of the 
planet. For Saturn we find a model cooling age ~20% longer than previous work. 
However, an additional energy source, such as that due to helium phase separation, is 
still clearly needed. For Neptune, unlike in work from the 1980s and 1990s, we match 
the measured T e s of the planet with a model that also matches the planet's current 
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gravity field constraints. This is predominantly due to advances in the high-pressure 
equation of state of water. This may indicate that the planet possesses no barriers 
to efficient convection in its deep interior. However, for Uranus, our models exacer- 
bate the well-known problem that Uranus is far cooler than calculations predict, which 
could imply strong barriers to interior convective cooling. The atmosphere grids are 
published here as tables, so that they may be used by the wider community. 

Subject headings: planets and satellites: atmospheres, interiors, individual: Jupiter, 
Neptune, Saturn, Uranus 



1. Introduction 



Planets cool as they age. For giant planets, what regulates the cooling of their mostly con- 
vective interiors is the radiative properties of the thin skin of atmosphere that rests atop the bulk 
of the planet's mass. The effect of atmosp h eres on giant pla net c ooling was first inve stigated in 
the mid 1970s, when braboske et all (fl975l) . lHubbardl Jl977h . and lPollack et all d 19771) computed 
the first thermal evolution models of Jupiter and Saturn that coupled the atmospheric and interior 
structures. 

The relationship between the atmosphere and interior has become even better appreciated 



since the discovery of hot Jupiters 15 years ago (IMayor & Oueloal 1995b . which orbit their parent 



stars at 5-10 stellar radii, and often intercept 10 4 - 10 5 more flux than Jupiter receives from the 
Sun . This incident flux drives the external radiative zone to press ures near 1 kbar in Gyr-old plan- 
ets dGuillot et all 19961 : ISudarskv et allEooilFortnev et al.ll2007h . which suppresses the transfer of 



riety of model atmosphere grids have been calculated for these close-in planets (e.g. 


Burrows et al. 


2003; 


Baraffe et al. 2003: Fortnev et al. 2007: 


Baraffe et al. 


2008; 


Ibgui & Burrows 


bood), which 


feature state-of-the-art treatments of chemistry (e.g. 


Lodders & Feglev 


2006 


), non-gray opacities 



(e.g. ISharp & Burrowsll2007l; iFreedman et al.ll2008l) . and radiative transfer 



These same atmosphere models have been honed ove r the past 10 years on excellent spectra 



of hundreds of brown dwarfs with T eff down to ~550 K (Allar d et al.ll200ll : iMarley et al. 



2002; 



Burrows et al.ll2006l : ISaumon et al.ll2006l : ISaumon & Marleyll2008l) . But these tools have not been 
applied to generate modern model atmosphere grids for use in evolution models of the solar sys- 
tem's giant planets. 

The first, and last, grid of non-gray radiative-convective model atmospheres to be tabulated 
in print for use in evolutionary calculations were presented in lGraboske et al] (|1975|) . These were 
computed for planetary effective temperatures of 20 to 1900 K, at only two surface gravities, that 
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of t he current Jupiter, a nd at 1/64 Jupiter's current gravity. More recently, iBurrows et al.1 (11997b 



and Baraffe etal 



( 2003 ) computed non-gray atmospher e models for use in the evolution of brown 



dwarf s and giant planets, including Jupiter and Saturn (|Hubbard et al 



20031). but the s e grid s are not available to the wider community. Both 
Burrows et al.1 ( 19971) did not treat irradiation from the Sun. 



1999 1: Fortnev & Hubbard 



Graboske et al 



Jl975l) 



and 



The purpose of the present paper is to provide radiative-convective model atmosphere grids 
for Jupiter, Saturn, Uranus, and Neptune, for the full ranges of surface gravities, effective tem- 
peratures, incident fluxes, and atmospheric compositions specific to these planets, as a function of 
planet age. This will greatly minimize any remaining uncertainty in cooling calculations related to 
the atmospheric boundary condition for these planets. Here we also apply these grids to investigate 
what effect they have on the cooling histories of the planets. 

We can briefly summar ize the main findings of evolution models of the major planets^. Go- 
ing back to iHubbardl (| 19771) . evolution models of Jupiter have yielded a planetary T e ff roughly 
cons istent with observations at an age of 4.6 Gyr, u nder the assumption of a fully a diabatic inte 



rior (Hubbard et al.lll999i : lFortney & Hubbardll2003l) . or with a small radiative zone (Gu illot et al. 
1995b. These models use an initial, post-formation T cfi - that is arbitrarily large (see e.g. lMarley et i 
2007bl) . on the order of 1000 K, implying that Jupiter was quite luminous at young ages, as is ex- 
pected from any model of its formation. Jupiter's evolution is considered fairly well understood. 

Fully convective homogeneous models of Saturn reach the planet's known T eS in only 2- 
3 Gyr, implying that Saturn is mu ch warmer than can be accommodated by the same kind of 
model that works well for Jupiter ([Pollack et al.1 1 19771) . The phase separation of helium from 
liquid metallic hydrogen, and its subsequent sinking to deeper layers (a differentiation process) 
has long been suspected of providin g Saturn's additional energy source (ISalpeterll 1973c IStevenson 



1975l : IStevenson & Salpete 



iHubbardetal.il 19991 : 



1977b 



the physics dMorales et al.ll2009TlLorenzen et al.ll2009l) and its implementation in cooling models 



a|). While this is wi dely assumed to be true, the exact details of 



Fortnev & Hubbard 



2003b has not yet been solved. 



Uranus and Nept une cooling mod els have their own problems , which are different than those 
of Jupiter and Saturn. IHubbardl (11978b andlHubbard & Macfarla nd (119800 computed the first ther- 
mal evolution models of these planets. iHubbard & Macfarlane could only find good agree- 
ment between model cooling ages and the age of the solar system, if the planets started out rel- 
atively cool, with post-formation T^s only ~ 50% greater than their current values. Put another 
way, the current planets were found to be underluminous compared to models that started with 
the arbitrarily large initial r eff values used for Jupiter and Saturn. The problem was especially 



'A nice revi ew of work until 1980, the early years of giant planet spacecraft observations and evolution modeling, 
can be found in Hubbard ( 1980) 
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pronounced for Uranus. 



Later work on cooling models can be found in lPodolak et al.l (| 199 11) and lHubbard et al.l (| 19951) . 
Again, both Uranus and Neptune appeared under -luminous compared to models, although the dis- 
crepancy was larger for Uranus than for Neptune. IPodolak et al.1 (1199 lb hypothesized that statically 
stable layers in the interior of the "ice giants," due to compositional layering, could lock in internal 
energy, leading to small intrinsic fluxes. They suggested that only the outer 40% and 60% in radius 
of Uranus and Neptune, respectively, were freely convecting. It is known th at the interior struc - 
tures of these planets do not ne ed to be partitioned into well-d efined layers (IMarley et al.N 19951) . 
More recently, and importantly, Fortney & Nettelmannl (120 lOh computed new cooling models of 
these planets. They used a new water equation of state (|French et al. I l2009h . and find that Neptune's 
current T eff can be matched, while the problem with Uranus remains, for models that feature the 
standard arbitrary hot start. The implication is that the argument for stable layers in Neptune may 
no longer be necessary. 

To summarize, current state-of-the-art models can reproduce the r eff of Jupiter and Neptune, 
while Saturn is over-luminous and Uranus is under-luminous. In this paper, in §2 we describe new 
model atmospheres for all of these planets, and apply them to new cooling models in §3. In §4 we 
give our conclusions and suggest future work. 



2. Model Atmospheres 
2.1. The model atmosphere code and validation 



An evolution calculation is dependent on a relation between th e intrinsic effective temperature 
Tint, surface gravity g, and specific entropy of the interior adiabat S (|Graboske et al.lll975l : lHubbard 
19771 : burrows et al.lll997|) . Here T int as defined as the r eff that the planet would have in the absence 
of solar insolation. With a third temperature, T eq , the r eff that the planet would have in the absence 
of internal energy, we have the relation r eff 4 =7j nt 4 +r eq 4 . Sometimes S is parameterized as T m 



the temperature of t he interior adiabat (which is isentropic) at P = 10 bars (IBurrows et al.lll997 : 
Hubbard et al.ll 1999T) . In practice, model atmospheres are calculated at a range of T mt and g points 
and the value of the r eff and T m are tabulated for each of these converged models. This collection 
of points in then inverted, such that for a given T 10 and g of a structural model, the r int and T eS are 
found by interpolation. While r eff can be most easily derived from observations, it is r int that is the 
parameterization of the flux from the planet's interior cooling. 

To create the grids we make use of a one-dimensional plane-parallel model atmosphere code 
that has been widely used for solar system planets, exoplanets, and brown dwarfs over the past 
two decades. The optical and thermal infrared radiative transfer solvers are described in detail 
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( Marlev & McKay 



2007; M arlev et al 



in iToon et al.1 (I1989I). Pa st applications of the model include Titan (IMcKav et al.lll989l). Uranus 



1999) 



19971: iMarlev et al. 



2007a; 



2002; 



Saumon & Marlev 



;as giant exoplanets ( Marleyl 199 8 ; Fortnev et al .12005 ; Fortnev & Marlev 



Fort nev et al. II2008L and brown dwarfs ( IMarlev et alJ 19961 ; burrows et al 



2008J). We use the correlated-k method for opacity 



tabulation ( boodyetalJI 1989b over 196 wavelength bins, from 0.26 to 325 /im for outgoing ther- 



mal radiation, and 160 wavelength bins, from 0.26 to 6 / im, fo r incident stellar radiation. Our 



extensive opacity database is described in lFreedman et al.l ( 2008b. We include t he opacity of neu 



tral atomic alkalis, which are prominent in b rown dwarfs ^Burrows et al.l 120001) . and also close a 



radiative region in Jupiter's deep atmosphere (IGuillot et al 



2004) 



We make use of tabulations of che mical mixing r a tios from equilibrium chemistry calcu 



lations of K. Lodders and col laborators ( Lodders 1999 ; Lodders & Fegley 20021; Visscher et al 



20061 ; lLodders & Fegleyl 120061) . We use the base protosolar abundances of lLoddersI (120031) . with 
Jupiter atmosphere models enhanced in heavy elements ("metals") by a factor of ~ 3 (specifi- 
cally [M/H]=+0.5), Saturn by a factor of 10 ([M/H]=+1.0), and Neptune and Uranus by a factor of 
~ 30 ([M/H]=+1.5). A 3 x enhancemen t in Jupiter is roughly consistent with observations by the 
Galileo Entry Probe (IWong et al.ll2004[). An en hancement of 10 x in Saturn is consistent with the 
CH 4 abundance deduced by LFlasar et al.1 (|2005l) via Cassini CIRS spectroscopy. An enhancement 
of 3 0-60 x in Uranus and Ne ptune is consistent with the methane abundance in these atmospheres 
(e.g. lGuillot& Gautiejl2009h . 



These are the first atmosphere grids for any of these planets to specifically include irradiation 
from the Sun. As mentioned, there h ave been only two prev ious non-gray atmosphere tabulations 
applied to evolutionary calculations, iGraboske et al.1 (| 19751) and iBurrows et al.1 (| 19971) . and only 
the former is publicly available. Both of these grids assumed solar abundances, no irradiation from 
the Sun, and were sparsely sampled in T sS and g. 

The utilization of no-irradiation atmosphere models rests on the assumption that one can 
assume that all flux absorbed by the plan et is absorbed into the planet's convective z one and directly 
adds to the energy budget of the planet (|Hubbardll 19771 ; LFortney & Hubbardll2003l) . Given that all 
the giant planets have stratospheres, due to absorption of incident flux in the radiative atmosphere, 
this assumption does not hold. In addition, one must assume a Bond albedo at every age, in order 
to determine the amount of absorbed incident flux. Typically, the current Bond albedo is used 
at all ages. Here we consistently solve for the deposition of stellar flux into the atmosphere as a 
function of wavelength and depth, and make no assumption regarding the Bond albedos. The 
at a given T mt and g is obtained based on the converged radiative-convective equilibrium model 
atmosphere temperature/pressure/opacity profiled As mentioned above, this improved method of 



2 The Bond albedo cannot be precisely determined from a ID model. A multi-D model must be used to investigate 
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treating the at mospheric boundary has been used for th e close-in "hot Jupiter" planets for several 
ye ars (see, e.g. iBaraffe et al.ll2003l : iBurrows et al.ll2003l for initial calculations) and was extended 
by Fortney et alT ( 2007 ) to models of EGP thermal evolution for planets from 0.02 to 10 AU. 



Here the derived P—T profiles are for planet-wide average conditions, with a Sunlight zenith 
angle of 60 degrees (/x= 0.5), and the incident flux cut by 1/2, due to the day /night difference 
(Ma rlev & McKavlll999l) . Using an analytical model for a semi-grey, non-scattering atmosphere 
(buillotlboiok we compared solutions obtained with that approximation to globally averaged solu- 
tions, with a greenhouse factor (i.e. the ratio of the infrared to the visible mean opacities) between 
10" 2 to 100. We found that mean deviations in the temperature profile between optical depths 
r = 10" 2 to 100 were smaller than 1.4%, with a standard deviation smaller than 3.2%, i.e. much 
smaller than other uncertainties. 

While the model atmospheres presented here are a significant advance over previous calcula- 
tions, some simplifications are still made. Most importantly, we neglect the opacity of condensate 
clouds and non-equilibrium hazes. For i nstance, the geome tric albedos of Jupiter and Saturn are 



depressed in the blue due to dark hazes (Karkoschka 
the photochemical destruction of methane ((Rages et aL 



1994) that are thought to be derived from 



1999). However, to our knowledge there is 



no adequate theory to predict how the optical depth and distribution of this haze may change as a 
function of solar luminosity, planetary T e s, and chemical mixing ratios, which change as the plan- 
ets cool. Even more importantly, equilibrium condensate clouds suc h as water or amm onia lead to 
large geometric albedos outside of methane absorption bands (e.g. ICahoy et al.ll2010|) . However, 
the cloud optical properties, particle sizes, vertical di stribution, and coverage fra c tion are diffi 



cult or imp ossible to predict from fir st principles. See lAckerman & Marleyl (120011) . ICooper et al. 



(|2003|) . and lHelling & Woitkd (120061) for modern cloud models for substellar objects. 



Of significant practical importance, and what validates our cloud-free assumption here, is that 
we reproduce the current temperature structure and r eff of the Jovian atmosphe re much better wit h 
cloud-free models than with models that include ammonia and water clouds (|Cahoy et al.ll2010l) . 
The cloudy models (which use the Ackerman & Marley cloud model) yield model atmospheres 
that are too reflective (a high Bond albedo) and are hence, too cool. This is because without dark 
photochemical hazes, the ammonia and water condensate clouds in the atmosphere model are more 
reflective than observed. Cloud layers in Uranus and Neptune reside deeper in the atmosphere, 
which makes this issue less important for their atmospheres. Looking back in time, when the 
planets will be warmer, points to a proper treatment of the time-varying Bond albedos being a very 
complicated problem. The cloud layers in Uranus and Neptune will be higher in the atmosphere. 
Jupiter and Saturn will lose their ammonia clouds, and possess a top layer of thick water clouds. 



the scattering of light over all phase angles (e.g. Marlev & McKavlll999 ; Cahov et"ap|2010 ) 
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There is much future work to be done in this area. 



Fortunately, for Jupiter and Saturn, the neglect of both the dark absorbing hazes, and rela- 
tively bright condensate clouds somewhat cancel out, and we are able to reproduce the temperature 
structure of the planets' tropospheres quite well. The exa mple of a cloud fr ee model of the current 
Jupiter, compared to data from the Galileo Entry Probe (ISeiff et al.lll998f) . is shown in Figured] 
The model T eS is correct to 0.2 K, and the 1-bar temperature (within the convective region) to 4 K. 
In Figure |2] we show a co mparison with the atmosphere of Neptune, which was probed via radio 
occultation by Voyager 2 (|Lindall 1 1 9920 . The comparison here is also favorable, with differences 
only on the order of several degrees. Since these planets cool very slowly at gigayear ages, a good 
match to T eii and 7i now also indicates a good match for the past several gigayears. 



2.2. The atmosphere grids 



Hubbard! (| 1 9771) provided an analytic fit to the original lGraboske et al.l(| 1975b non-gray (mean- 
in g frequency-dependent opacities are used) grid of m odel atmosphere calculatio ns. As discussed 



in 



Hubbard & Macfarland (I19801) . iGuillot et all (11995b . and lSaumon et all (11996b . this grid can be 



described by a function of the form 



V -0.167^1.243 



(1) 



where T is the temperature at a reference pressure Pp within the conve ctive region of the atmo- 
sphere, such as 1 or 10 bars, and K is a constant. ISaumon et al.l (| 19961) find that T =\0 bar and 
K=3.36 provides a reas o nable fit below r eff of 200 K. For Uranus and Neptune evolution models, 
Hubbard & Macfarland (|1980|) . in order to better match the planet's atmospheres, suppressed the 
gravity dependence entirely and set K=l, with To = 0.750 bar for Uranus, and P = 0.861. So, in 
general, a variety of fits, some based on the original calculations, some based on ad-hoc modifica- 
tions, are available. Clearly new work in this area is warranted. 



The lBurrows et al.l (|1997b grid, slightly updated in lHubbard et al.l (|1999|) . spans a very wide- 
range of T eS , g, and Tio, with a mix of gra y atmosphere s at hig h T eS , and non-gray atmospheres 
at low r eff , suitable for giant planets. The iBaraffe et al.l (|2003b grid, fully non-gray, is similarly 
expansive, to treat the lowest mass stars, brown dwarfs, and planets. 

We have computed pressure-temperature (P-T) profiles for Jupiter and Saturn at ~50 T eS /g 
pairs, enough for detailed coverage of the evolution of Tj nt , g, and S, over their evolution. For 
Jupiter, ^ ranges from 89 K to 1200 K, while the gravity range covers 0.12 to 1.1 x Jupiter's 
gravity. For Saturn, T int ranges from 69 K to 650 K, while the gravity range covers 0.12 to 1.1 
x Saturn's gravity. For both planets, these calculations were done at the current value of the 
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incident solar flux, a s well as 0.7 of this value so that the evolution of the solar luminosity can 



be incorporated (e.g. iHubbard et al.lll999|) . These are included as Tables 1 and 2, for Jupiter and 
Saturn, respectively. The grids for 1.0 L & for Jupiter and Saturn are shown in Figures [3] and HI 
respectively. They show T mt and T efi on the x-axis, and T 10 on the y-axis, for 8 different surface 
gravities. Low gravities cover high 7j nt (young ages) and high gravities cover low r int (old ages). 
An inset shows T e s vs. T w for the current and lower solar luminosity, at very low 7^. When 7i nt is 
large, then the planet's energy budget is dominated by its own intrinsic luminosity, and r eff is only 
negligibly larger than 7j nt . However, at lower values of T m , eventually there is a clear separation 
between 7} nt and T eS , as expected. 

For Uranus and Neptune, we have computed four separate grids, at four different values of 
the incident flux. The 7j nt ranges from 27 to 217 K, while the gravity range covers 0.45 x Uranus' 
gravity to 1.1 x Neptune's gravity. Two values of the incident flux are for the current values for 
Uranus and Neptune, at 19.2 and 30.1 AU from the current Sun, respectively. The other two are 
at 1.8 x the current flux received by Uranus, and 0.12x the current flux received by Neptune. 
This very wide range allows for the inclusion of a time variable solar luminosity a s well as the 



possi bility of dramatic changes in the orbital distances of the planets with time (e.g. iGomes et al. 



20051) . The grids at Uranus' and Neptune's current flux levels are compiled in Table 3 and shown 
in Figures [5] and [6l respectively. At high r eff the two grids appear quite similar, but for Uranus' 
larger incident flux, the wider division between 7i nt and T e & at low 7] nt , is readily apparent. For 
nearly all models of the four planets, the atmosphere becomes convective and stays convective by 
10 bars. For a handful of models for Jupiter, the deep convective adiabat was not reached until 
~ 15 bars, but the tabulated valu e of T w is modified to be consistent with the entropy of the deep 



adiabat (see lBurrows et al.lll997|) . 



The model atmosphere grid, and its relation between S, 7] nt and r eff , comes in through the 
energy conservation equation 

dL _ ^dS ^ 
dm dt 

where L is the luminosity, m is the mass variable, T is the temperature in a mass shell, S is the 
specific entropy of that shell, and t is the time. After explicitly defining the dimensionless mass 
shell variable m as 

If 2 

m= — / Aixr' dr'p(r'), (3) 
M Jo 

we can rewrite Eq. © in terms of the time step dt, as 



M 



dt = -— TdSdm, (4) 
L Jo 

where dt is the time step between two models that differ in entropy by OS. The luminosity L 
extracted from the planet is 4nR 2 aT^ v where the value of T{ nt at a given S is given by the model 
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atmospheres. While at young ages 7i nt ~ T e s, at old ages, when T[ nt can be appreciably smaller than 
T ef f, interior cooling can be quite slow. After presenting our new cooling calculations in the next 
section, in §13.31 we will examine the affects of the new atmosphere grids, compared to previous 
tabulations. 



3. Cooling Calculations 
3.1. Jupiter and Saturn 



Using these new model atmosphere grids, we can calculate the cooling history of the major 
planets. Cooling calculations f or these planets have b een published by many authors. Recent work 
for Jupiter and Saturn includes iHubbard et al.1 dl999|) . who explicitly showed the effects of a faint 
young Sun, and who included some limiting cases of additional interior en ergy due to the phase 
separa ti on of helium from l i quid m etallic hydrogen at megabar pressures ([Stevenson & Salpeter 
1977a|) . iFortney & Hubbard! (I2003|) looked at the evolution of Saturn and tested a num ber of previ- 
ously published phase diagrams for H/He phase separation. ISaumon & Guillotl (120041) investigated 
the sensitivity of Jupiter cooling models to various hydrogen equations of state, which can pre- 
dict quite different te mperatures in the deep interior of the planet. All of these models used the 
Burrows et al.l (| 1997f) atmosphere grid. 



In Figures[7]and[8]we compare Jupiter and Saturn evolution models to those of lFortney & Hubbard 



(120031) . The old and new models make the same assumptions regarding solar luminosity and plan- 
etary structure, so the only differences are due the model atmospheres grids. For both planets, a 
helium mass fraction (Y), relative to hydrogen, of 0.27 is used in the adiabatic H/He envelope. 
An rocky (zero-temperature ANEOS olivine) core mass of 10 M e is used for Jupiter, and 21 
M ffi for Saturn. Within the H/He envelope, the zero-temperature ANEOS water EOS is used to 
mix in a mass fraction of 0.059 and 0.030 of "metals" in Jupiter and Saturn, respectively. In 
Fortney & Hubbard! (|2003l) these choices reproduced the current radius and axial moment of iner- 
tia (obtained fro m more detailed structure m odels) at the time the planets cooled to their known 
r e ff values. As in Fortnev & Hubbard! (120031) the heat c ontent of the rock and water are neglected 



in the evolution calculation. IFortney & Hubbard! (120031) use a Bond albedo 0.343 was assumed for 
both planets at all ages, while now we use the self-consistent grids. The cooling times are modestly 
prolonged, as will be discussed in detail in § 13.31 

A particular interesting difference shown in Figure [8] is the larger radii for the new planet 
models at young ages, which is most pronounced in Saturn. This is a manifestation of the arbitrary 
initial condition (a hot, ~ 3 Rj, adiabatic sphere) along with the slowed cooling in the 700-400 K 
7j nt range compared to previous models (See §13.31) . The initial conditions for cooling are tied to 
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details of the energetics of planet formation, and are not well understood (|Marley et al.ll2007bl) 



The revised cooling age for Saturn does little to change the l ong-standing proble m that the 



planet is much more luminous than homogeneous models predict (|Pollack et al 



1977 ). This has 



long been attributed to the phase separation ("demixing") of helium from liquid metallic hydro 
gen. This helium solubility is thought to be minimized at pressures of several megaba rs, at pres 
sures where the (gradual?) dissociation and ionization of hydrogen is nearly completed (|Stevenson 
1975b . The phase diagram of H/He mixtures is beyond the realm of curre nt experiment, but re 



cent advances in fi rst principles calculations of the H/ He phase diagram (ILorenzen et al. 



Morales et al. Il2009h should be tested in evolution models (IHubbard et al. 1 19991 : 



200 



9: 



Fortnev & Hubbard 



2003|) . to see if the additional energy source from this "heli um rain" can explain Saturn 's luminos- 
ity. Previously calculated phase diagrams were tested in iFortney & Hubbard! (|2003|) . and were 
found to not explain Saturn's thermal evolution. A complication that must be addressed in the fu- 
ture is whether the deep interior temp erature gradient becomes dramatically super-adiabatic in the 
face of helium composition gradients ([Stevenson & Salpetenll977al : LFortney & Hubbardll2003l) . 



The model calculations for Jupiter, which yield an age of 5.3 Gyr, rather than 4.6 Gyr, could 
have important implications for the p lanet. Jupiter's atmo sphere is clearly depleted in helium, 
according to Galileo Entry Probe data (|von Zahn et al.ll 19981) . which is a strong indication that he- 
lium phase separation has occurred in this planet. Furthermore the atmosphere's depletion in neon 
(|Mahaffy et al] 120001) is strongly suggestive of h elium demixing as well, as neon is expected to 
prefer entially dissolve into helium-rich droplets (|Roulston & Stevenson! 1 1995t IWilson & Militzei 
2010j). The inclu sion of the additiona l energy release due to this demixing will further prolong 
Jupiter's cooling (IHubbard et al.ll 19991) . leading to a worse match with observations. 



However. IS aumon & Guillotl (|2004l) have investigated cooling models of Jupiter with a variety 
of hydrogen EOSs, that predict a wide range of temperatures in Jupiter's deep interior, which led 
to evolutionary ages for homogeneous m odels as short as 3 G yr. If the deep interior temperatures 
in Jupiter are lower than found with the ISaumon et al.l (I1995|) EOS used here, then it is possible 
that a combination of the model atmospheres presented here, helium rain (which prolongs the 
evolution), and colder interior temperatures (which quickens the evolution) c ould yield a good 



match to observations. R ecent work on the hy drogen EOS, both theoretically (INettelmann et al 



20081 : Militzer et al.ll2008l) . and expe rimentally (IHolmes et all 1 19951 : ICollins et al 



temperatures lower than predicted by lSaumon et al.l(|1995l) . so this avenue is plausible. 



20011) . do yield 
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3.2. Uranus and Neptune 



The long-standing problem for Uranus and Neptune has been that both of these planets 
are colder at the present day than cooling models pr edict. This is a reve rse o f the situation 
for Saturn. This issue is discussed in some detail in 



Podolak et al. 



(1199 II) and Hubbard et al 



(11995b . and in the general literature in iHubbard & Macfarland (I1980I) . In order to understand 
how advances in input physics over the past 30 years in the EOSs affect the thermal evolu- 
tion of these planets, we hav e computed a set of evolution models that use the physics of the 
Hubbard & M acfarlane ( 198 01) models, and compared them to our new calculations. The models 
presented in Hubbard & Macfarlane d 1980) hav e three distinct, adiabatic, layers. The H/He enve- 
lope uses the E OS of battery & Hubbard! d 1976), the "icy" layer uses a mixture the H 2 0, NH 3 , and 
CH 4 EOS from lzharkov & Trubitsvnl d 1978k and the rocky core ( a mixture of silicon, magnesium , 
iron, oxygen, and sulfur) is also from lzharkov & Trubitsynl (| 19781) . IHubbard & Macfarland (119801) 
also make estimates for the specific h eat capacity of the icy and rocky layers. Our implementation 
of the IHubbard & Macfarland (|1980h models agree well with the original, particularly for Uranus, 
and yield slightly shorter cooling times for both planets (dashed curves in Fig. |9]), compared to 
their work. 

To investigate how updated EOSs affect the evolution, we can cr eate cooling models that 
use the same ice-to-ro ck ratio as used in IHubbard & Macfarland (|l980|) . 2.71-to-l, and the same 
Graboske et al] (| 19751) atmosphere grid. But instead we use updated EOSs for all three layers. 



These are the 



Saumon et al. 



(1995) EOS for H/He and the Sesame EOSs of "water 7154" and 



"dry sand" ( Lyon & Johnson 19921) for water and rock, respectively. These tabulated ice and rock 
EOSs include calculations of the density- and temperature-dependent free energy at every temper- 
ature/density point, so no assumption must be made for the average heat capacity. In Figure [9] we 
compare cooling tracks with old and new EOSs. Following previous work, we plot these tracks 
backwards in time from the current day, to see what initial values of T e s could explain the current 
planets. The updated EOS yield much faster evolution for both planets (dotted curves). The change 
for Neptune is enough to allow the planet have an initially "hot start." However, Uranus models 
must start at a very low r eff to explain the planet's current low T eff . Using these same new interior 
models, but with our new model atmosphere grids, yield the thick solid curves in Figure |9j These 
lead to slower cooling for both planets, most dramatically for Uranus. Neptune's evolution is still 
approximately co nsistent with an arbitrari l y hot start at formation. This last finding agrees with the 
work published in lFortney & Nettelmannl (|20 100 . 



While the 3 -layer models presented in IHubbard & Macfarland (|l980l) were at the time con- 
sist ent with observational c o nstrai nts on their interior structure, that is no longer the case. Neither 
the IHubbard & Macfarland (|l980l) models, nor our new implementation of their 3 -layer models 
(with the revised EOSs), are consistent with the gravity fields of these planets. To further expand 
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our treatme nt of Uranus and Neptu ne, we calculated new structure m odels, using the methods 
described in lNettelmann et al.1 (I2008|) and lFortney & Nettelmannl (|2010|) . We then investigated the 
thermal evolution of these models that are also consistent with the constraints on current structure^ 



Like in iFortney & Nettelmannl (|2010l) . these models also use a three-layer structure, but in- 
clude some water mixed into the H/He layer, and some H/He mixed into the water layer, above 
the rock core. The outer layer is predominantly H/He with a helium mass fraction, relative to hy- 
drogen (Y), of 0.27 (the H-REOS and He- REOS equations of state iNettelmann et al.lliooih with 
a 0.30 mass-fraction of water (the EOS of French et al.l 120091) . The middle layer is mostly fluid 
water (mass fraction 0.878, beginning at 0.20 Mbar for Uranus, and 0.852, beginning at 0.10 Mbar 
for N eptune), with a small adm ixture of H/He. Both of these layers are adiabatic. The core is 
rock (|Hubbard & Marleyll 19891) with a mass of 1.51 M e for Uranus and 2.85 M m for Neptun e. For 
the evolution models, the rock core uses a radioactive luminosity from lGuillot et al.l (I1995|) and a 
specific heat capacit y of 1 Jg-'K" 1 . The ful l allow ed range of U ranus/Neptune interior composi- 
tions are explored in IFortney & Nettelmannl (120101) . The fit to the lGraboske et al.1 (119751) grids use 
Ti=73 K (AT=1.418) and Ti=78 K (^=1 .571) for Uranus and Neptune, respectively, including the 
gravity dependence. 

Evolutionary calculations for these model planets are shown in Figure [T0l We can repro- 
duce well the current T e $ of Neptune (blue models) with our new model atmosphere grid {solid 
cu rves), as well as the o l der gr id {dashed). The latter models agree well with work published 



in 



Fortney & Nettelmannl (|2010h which shows that the cooling times are insensitive to details of 



structure assumptions within the range allowed by the observational uncertainties. Therefore, there 
is a plausible consistency for the planet: the current interior structure and cooling history can be 
matched by one model. As shown in Figure [HE the mismatch for Uranus becomes appreciably 
larger with the new model atmospheres, in agreement with Figure [9] For Uranus in particular the 
model indicates a very slow change in r eff with time in the current era, due to the larger incident 
flux and higher T eq than for Neptune. Therefore small changes in the model atmosphere can lead to 
dramatic changes in cooling times, as seen in Figures [9] and [T0l Another manifestation of Uranus's 
slow cooling is shown by the thin red curve in Figure[9| T his shows that a tiny cha nge in the current 
ri nt , to the maximum allowed by observation {la error bar. lPearl & Conrathlll99ll) can dramatically 
alter the calculation of the past cooling history. 



These particular curves use T\ = 73 and 78 K for the current Uranus and Neptune, to better allow the static 
structure model to match the T e g at the current time, at the expense of a slightly worse match to the observed T\. 
However, changes in T\ by a few degrees in either direction at the current time, which adjusts the interior structure, 
has little effect on the cooling history. 
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3.3. Effects of the Atmosphere Grids 



As shown in Figures ITlfTOl a general finding of our work is that the cooling of our solar 
system's giant planets is slowed with the new model atmospheres. The main reason is higher 
atmospheric opacity, due to improvements in opacity datbases, especially at high temperature, as 
well as the inclusion of opacity sources not previously known in 1975 or 1997. Here we investigate 
the reasons for these differences in atmospheric structure and cooling. We will first examine Jupiter 
and Saturn. 



It is important to note that the agreement between the iBurrows et al.1 dl997|) grid and our 
grids are best at current (low) 7j nt value s. This is not necessarily surprising. Both works use 
the same model atmosphere code (that of lMarley et al.lll996l : lMarley & McKaylll999l) . However, 
the opacity databases used have changed significantly over the past 14 years, most importantly 
at high Tim values. At low temperature, the only remaining opacity sources are methane vapor 
and H 2 collision induced absorption. At warmer temperatures, observations of numerous T-type 
brown dwarfs have necessitated dramatic revisions to model atmospheres s ince 1997, including 
the inclusion of i mportant new opacity sources, such as the alkali metals (IBurrows et al.l 12000 : 
Allard et ap200 1 ) . In addition, high-temperature molecular opacity databases, which were in their 
infancy in 1997, are now becoming mo r e com plete, which generally leads to higher opacities. 
Compared to Figure 2 of IBurrows et al.l (| 19971) we find much smaller detached radiative zones 
below the photosphere, over a narrower range of 7j nt . (In fact, for the Jupiter grid specifically, only 
the 450 K and 600 K models possess them.) 

At contant T mt , a smaller or nonexistent detached radiative zone leads directly to a lower 
specific entropy adiabat, as shown in Figure [TT]for models at 596 and 417 K, representing Jupiter 
at ages of 10 and 32 Myr, respectively. In black we show the standard 3 x solar metallicity models. 
In orange we show the same calculation, but with 1 x solar metallicity. The photospheric pressure 
(where T=T e s) moves to higher pressure, due to the lower opacity, yielding a colder (lower specific 
entropy) adiabat. The extent of the detached radiative zone is only marginally affected — the general 
trend towards smaller radiative zones at lower r eff is not disturbed. In magenta we plot the same 
1 x solar metallicity models, but with the Na and K alkali opacity removed, to show their affect. 
The radiative zones are larger in vertical extent, with a shallower temperature gradient. This shows 
that alkali opacity is a strong cont r ibutor towards closing the radiative zone, in a manner similar 
to that suggested by iGuillot et al.l (|2004l) for Jupiter's current putative radiative zone at similar 
temperatures (1000-2000 K). However, we note that even with alkali opa cities removed, we ar e 
still unable to match the large extent of the radiative z ones from Figure 2 of IBurrows et al.l (| 19971) . 
which shows that overall molecular opacity updates ( Freedman et al. 2008|) since that time also 
contribute to the smaller radiative zones found today. 



We now turn to the affect on evolution. Given our discussion above, a given drop in S between 
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two identical structure models will occur at a lower T[ nt for a fully convective atmosphere, which 
leads to a longer time step in Equation [4] since this step goes at T int to the fourth power. Another 
way to look at cooling differences between the old and new atmosphere grids can be gleaned from 
Figure [121 We plot planet radius vs. T inU as the radius can serve a proxy both for interior specific 
entropy (parameterized by T lQ ), and we can visually examine the surface gravity and 7j nt changes 
experienced in the model. On this plot regions where the new atmosphere grid lead to slower 
cooling are where the slope dR/dT inU our proxy for dS/dT mt , is steeper for the black curves than 
that of the red curves. For a give change in r int , they have a larger AS, meaning a longer time 
step from Equation[4] For example, a comparison with 7j nt vs. time in Figure [7J shows that the new 
models lead to slower cooling from T int of ~800 K to the present 7j nt , although the differences are 
extremely small at 7] nt below 200 K. 

In Figure [13] we show model P—T profiles for Jupiter's atmosphere, corresponding to the 
new cooling curve for Jupiter in Figure [7J While the early evolution is highly uncertain (e.g. 



Marley et al.ll2007bl) . we find that at 1 Myr Jupiter's atmosphere had only one (deep) convective 



zone, while a second, detached convective zone appears for a few tens of millions of years, from 
perhaps ~5-50 Myr. We also find that Jupiter's water clouds formed at an age of ~ 30 Myr, while 
the ammonia clouds began to condense at ~ 1 Gyr.. 

Figure [14] shows that the Uranus and Neptune evolution models generally show an effect 
similar to that seen for Jupiter and Saturn. Here the new grids always lead to a smaller 7] nt at 
a given radius. Hence the cooling is always slower when using the grids, compared to those of 



Graboske et al.l (|1975l) . It is clear from Figure [14] that the apparently small difference between the 
grids at low 7j nt becomes magnified because the cooling is so slow. Even at higher 7i nt , farther 
back into the past, the differences betw een the two grids is lar ger than in Neptune. There, the good 
agreement at low 7] nt is an effect of the lGraboske et al.l (| 19751) grids being ad-hoc "tuned" to agree 
well with the current Uranus and Neptune, with P r eference pressure being chang ed from 1 bar to 
0.750 bar for Uranus and 0.861 bar for Neptune (see lHubbard & MacfarlaneHl980|) . As one moves 
to higher T[ nt , there is no reason to expect this tuning of the grid to hold, so the old and new grids 
diverge. 



4. Discussion and Conclusions 



Our own giant planets serve as our calibrators for the evolutionary theory used to under- 



stand the thermal evolution of extrasolar giant planets (e.g. iHubbard et al.l 120021 : iGuillotl 12005 



Fortney & Nettelmannll2010l : iBaraffe et al.ll2010l) . As new equations of state come online, better 
understanding of interior energy transport develops, or new calculations of atmospheric structure 
or opacities are made, the problem of the cooling of the giant planets needs to be revisited. We 
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have investigated "standard" cooling models for these planets, where the planetary luminosity is 
predominantly due to the slow release of remnant formation energy, within adiabatic layers beneath 
the radiative atmosphere. 

If our calculations are taken at face value, we can robustly conclude, as have authors before 
us, that Saturn i s quite overluminous, necess i tating a large role for heliu m phase separation over 
the past 2 Gyr (IStevenson & Salpeterl 1 1977al : LFortney & Hubbard! 120031) . For Jupiter, the model 
cooling age is now modestly too long, which may suggest that the inte rior temperatures are c older 
than yielded by our models, as has been hinted at previously. Although lHubbard et al.1 (119991) have 
shown that cooling ages for Jupiter of 3.5 Gyr are possible in the absence of incident solar flux 
(a Bond albedo of 1) for this H/He EOS, extreme atmospheric reflectivity does not seem to be a 
realistic pathway to faster cooling. 



Like others, we find that Uranus is quite underluminous. IPodolak et al.l (| 199 II) have suggested 
that a statically stable interior, unable to convect due to composition gradients, may be the best 
explanation. However, these same models for Neptune can match the planet's T e f{ at age of 4.5 
Gyr. The improvement is based predominantly on modern updates to the EOSs of constituent 
materials. It appears that whatever deep interior complications arise in Uranus do not arise in 
Neptune, or to a much smaller degree. This could have profound implications for convection in 
these planets, and in pa rticular the generation of the se planets' non-axial non-dipolar-dominated 
magnetic fields (see, e.g lStanley & Glatzmaienl2009[ for a review of dynamo modeling). 



Us ing the interior statically stable geo metry suggested by lPodolak et al.l(|l99ll) and lHubbard et al 



(11995b . IStanley & Bloxhaml (|2004l |2006|) were able to reproduce the main features of the mag- 
netic fields of Uranus and Neptune by hypothesizing dynamo action in only an outer shell in both 
planets. However, we find that such a picture may not be viable for Neptune. More recently 
Gomez-Perez & Heimpell (|2007|) have constructed 3D dynamo models with radial variable elec- 
trical conductivity, and found similar good agreement with observations without adhering to the 
statically stable interior geometry. More work in this area is certainly needed. 

The summary figure of the models presented here — the luminosity of the giant planets with 
age — is shown in Figure [151 Of course the luminosities at young ages for all planets a re uncer- 
tain, and depend strongly on the details of the formation process (|Marley et al.ll2007b|) . Lower 
post-formation luminosities are certainly realistic, although at gigayear ages the details are not 
important. In addition the cooling curve of Uranus should be regarded with suspicion. 

One could certainly envision more complex models of the interiors of these planets, particu- 
larly for Uranus and Neptune. We have modeled the "fluid ice" component of these planets with 
a water EOS, but some previous authors have chosen instead to use mixtures of water, ammonia, 
and methane. Shock experiments have studied a C-N-O-H mixture called "synthetic Uranus" for 
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use in ice giant modeling (INellis et al.ll 19971 : iHubbard et al.lll99l|) . The relative amounts of C, N, 
and O in these planet is not constrained by structure models or formation th eory. Also, one mus t 
remember that the high-pressure EOS of planetary material are uncertain. IBaraffe et al.1 (I2008|) 
have quantitatively explored the use of different EOS for water, rock, and iron on the structure and 
evolution of Jupiter- and Neptune-class exoplanets. The evolution of Uranus and Neptune could 
be most affected by EOS uncertainties. 

We also caution that the good agreement with observations for the Neptune cooling models 
could be a coincidence. Theory and experiment have probed the phase diagram of pure carbon, and 
have shown that it is solid diamond at Neptune-interior P - T c onditions. A rain of solid diamond 
has been suggeste d for the interiors of Uranus and Neptune (|Rossl Il98ll : iBenedetti et al.1 Il999i : 
Eggert et al] 1201 01) . which could, at least in principle, be an additional energy source that prefer- 
entially powers Neptune more strongly than Uranus, to explain their dichotomy. These avenues 
should be explored in the future. 

There are a few paths towards improving the atmosphere grids presented here. For Jupiter and 
Saturn, one could investigate the reduced He/H ratio as the planets' age, which would affect the 
hydrogen collision induced absorption (CIA) that is an important infrared opacity source in these 
atmospheres. We recom mend updating CIA opacity in general, as the state-of-the-art in H-CIA 
opacity calcul ations (e.g. iBorvsow |2002|) are now showing some mismatches in modeling brown 
dwarf spectra (|Cushing et al.ll2008|) . The inclusion of condensate clouds, either from equilibrium 
chemistry, such as ammonia and water, or a methane-derived photochemical haze, are important in 
accurately modeling the energy balance and temperature structure of these atmospheres. Whether 
this could be understood well enough to predict with confidence the effects of clouds over the 
range of past 7^, atmospheric chemical mixing ratios, and surface gravity, in a ID planet-wide 
average model atmosphere, is an open question. It seems likely that further improvement towards 
understanding the cooling of the planets will likely come from work on the EOS of planetary 
materials. 
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T (K) 

Fig. 1. — Comparison of a model of Jupiter's current equatorial P-T profile (dotted black line) to 
that obtained for Jupiter by the Galileo Entry Probe (thick gray lines). The model parameters are 
shown on the plot. The value of T int is chosen to match the value of the intrinsic flux measured by 
Voyager 2. The resulting T e g is 124.6 K, the observed value within the observational error bars. 
The model planet's 1-bar temperature is cooler by 4 K than that measured by the probe. 
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Fig. 2. — Comparison of a model of Neptune's current P-T profile (dotted black line) to that ob- 
tained for Neptune via radio occultation by Voyager 2 (thick gray line). The model parameters are 
shown on the plot. The value of T{ nt is chosen to match the value of the intrinsic flux measured by 
Voyager 2. The resulting 7^ is 59.2 K, the observed value within the observational error bars. The 
model planet's 1-bar temperature is warmer by 5 K than that measured by Voyager 2. Extrapolating 
to the 10-bar level, the temperature appears to be a good match. 
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Fig. 3. — The Jupiter model atmosphere grid at 1.0 L & . Plotted on the x-axis is the T eff (solid 
lines) and 7^, (dotted lines) vs. T\q. Lower gravity models (yellowish-green, orange) cover only 
high T\q and T e ff. Higher gravity models (red, bright green) cover only low T l0 and r e ff. The effects 
of irradiation from the Sun (seen as a higher r eff than r int ) are only apparent at low temperatures. 
Inset: T eS vs. T w for the 1.0 L Q (solid,with squares) and 0.7 L (dash-dot, with triangles) grids. 
This shows that the effects of the lower solar luminosity do have small, but quantitatively important 
effects at low temperature. 
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Fig. 4. — The Saturn model atmosphere grid at 1.0 L & . Plotted on the x-axis is the r eff (solid lines) 
and Tim (dotted lines) vs. T w . Lower gravity models (yellowish-green, orange) cover only high 
7io and T eff . Higher gravity models (bright green, red) cover only low 7\o and T eS . The effects 
of irradiation from the Sun (seen as a higher r eff than r int ) are only apparent at low temperatures. 
Inset: T eS vs. T w for the 1.0 L (solid, with squares) and 0.7 L Q (dash-dot, with triangles) grids. 
This shows that the effects of the lower solar luminosity do have small, but quantitatively important 
effects at low temperature. 
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Fig. 5. — The Uranus model atmosphere grid at 1.0 L Q . Red is the lowest gravity, 4.0 m s" 2 , while 
yellow is the highest, 12.6 m s~ 2 . The x- and y-axes are the same as in Figures [3] and HI meaning 
the r e ffS are solid lines and 7i nt s are dotted lines. 
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Fig. 6. — The Neptune model atmosphere grid at 1.0 L Q . Red is the lowest gravity, 4.0 m s" 2 , while 
yellow is the highest, 12.6 m s" 2 . The grid is generally similar to that for Uranus, shown in Figure 
|5l but since the incident flux upon the planet is smaller, one generally find a smaller T e s for a given 
7] nt . The r eff s are solid lines and 7] nt s are dotted lines. 
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Fig. 7. — Thermal evolution models of Jupiter and Saturn. The upper panel shows the full evolu- 
tionary history, while the lower panel more clearly focuses on the current era. Solid curves are T e s, 
while dotted curves are 7] nt . The thick black curves use the model atmosphere grids presented in 
this work, including the time- variable luminosity of the Sun, a linear increase from 0.72 L Q at time 
zero. The thinner red curve s utilize t he same interior structure and solar luminosity, but instead 
use the lBurrows et alJ (| 19971) grid (see LFortney & Hubbardll2003r ). The upper curves are for Jupiter 
while the lower curves are for Saturn. The current uncertainty is the T eff of each planet is repre- 
sented by the size of the square at 4.56 Gyr in the lower panel. The Jupiter model is over-luminous 
by 10%, with a mean radius of 69,400 km. The evolution of Saturn does not quite reach the current 
time, but a slight extrapolation off of the atmosphere grid shows the Saturn model has only 63% of 
the actual planet's luminosity, with a mean radius 56,200 km. 
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Fig. 8. — Radius evolution of the models of Jupiter and Saturn shown in Figure [7] . The black 
curves use the new model atmosphere grids presented in this w ork. The red curves utilize the 
same interior structure an d solar luminosity, but instead use the iBurrows et all (|1997h grid (see 



Fortnev & Hubbardll2003 ). 
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Fig. 9. — Thermal evolution models of Uranus (red) and Neptune (blue) for models that feature 
three distinct homogeneous layers of H/He, fluid ices, and rock. The models are constructed to 
match the of both planets at the present day, and the evolution is followed backwards in time. 
The thick grey line indicates the form ation of the solar system. The dashed cur ves are cooling 



calcul ations performed in the manner of lHubbard & Macfarlanel (|1980f) . and use the lGraboske et al 



(|1975|) model atmosphere grids. The dotted curves use the same model atmospheres, but updated 
water and rock EOSs, which dramatically shorten the cooling times for both planets. The thick 
solid curves use the same updated EOSs as the dotted-curve models, but also use our new model 
atmospheres. The thin solid red curve, for Uranus, changes the thick- solid-red model only by 
adjusting the current T im to the upper boundary of the la error bar. 
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Fig. 10. — Thermal evolution of Uranus (red) and Neptune (blue) for models that match the current 
gravity field constraints. The thick solid curves use the mo del atmosphere grids presente d in this 
work. The dash ed curves are similar t o those presented in iFortney & Nettelmannl (|2010h . which 
used a fit to the iGraboske et al.1 (| 19751) grid. When the T eS of each planet reaches the observed 
values, the interior structure of each planet is consistent with the gravity field. At 4.56 Gyr, the 
radii of Uranus and Neptune are 25730 km and 24670 km, respectively. Inset: A zoom-in on the 
current era, which shows that the T e s of Neptune is well- reproduced, bu t the model Uranus is far 
too luminous. Measured T eS and error bars are taken from lGuilloti (I2005|) . 




T (K) 

Fig. 11. — Model pressure-temperature profiles of Jupiter at 417 K (g = 20.1 m s" 2 ) and 596 
K (g = 17.7 m s" 2 ). Thick lines indicate convective energy transport, while thin lines indiciate 
radiative energy transport. Cloud condensation curves are in dotted black. The equal-abundance 
curves of CH 4 /CO and NH 3 /N 2 are dashed. The solid black curves represent the atmosphere of the 
Jupiter model at 32 and 10 Myr, and feature 3x solar metallicity. The orange curves are the same, 
but feature solar metallicity. The magenta curves are also solar metallicity, but with the opacity 
of neutral atomic Na and K atoms removed. From black to orange, decreased metallicity leads to 
decreased opacity, a higher pressure photosphere, and lower specific entropy adiabat. The removal 
of the alkali opacity leads to a detached radiative zone that is larger in extent with a shallower 
temperature gradient, again leading to a cooler adiabat. 
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Fig. 12. — Radii of models of Jupiter and Saturn from Figure[7]as a function of 7j nt . (1 Rj is 71492 
km, Jupiter's equatorial radius.) Recall that radii, S, and T w are all monotonic functions of each 
other. As the planets cool, their radii contract, and S and T w decrease. The open circles are time 
steps at ages of 10 6 , 10 7 , 10 8 , 10 9 , and 10 10 years. Steeper dR/dT int slopes lead to slower cooling. 
(See text.) For instance, a comparison with Figure [7] of 7; nt for Jupiter shows that the new (black) 
models, compared to the red models, cool more quickly at T int > 800 K, but more slowly below 
800 K. Below 200 K, the differences are quite small. The same general picture is true for Saturn. 
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Fig. 13. — Model pressure-temperature profiles of Jupiter from log t=6 to 9.5 years, in steps of 0.5. 
The thin blue curves are radiative regions while thick black curves are convective regions. Cloud 
condensation curves are in dotted black. The equal-abundance curves of CH 4 /CO and NH 3 /N 2 are 
dashed. The rightmost profile is T eS = 976 K, g = 1 1 .4 m s~ 2 at 1 Myr. The profiles 3rd and 4th from 
the right (596 K, 17.7 m s" 2 , 10 Myr, and 417 K, 20.1 m s" 2 , 32 Myr) have a detached convective 
zone. From the models, it appears that Jupiter's water clouds formed at an age of ~ 30 Myr, while 
the ammonia clouds began to condense at ~ 1 Gyr. 
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Fig. 14. — Radii of models of Uranus and Neptune from Figure |9] as a function of T int . Recall 
that radii, S, and 7\ all monotonic functions of each other. Like Figure [H this plot is viewed as 
going backwards in time. As the planets are warmer in the past, their radii enlarge, and S and T w 
increase. The closed circles are time steps, equally space in 1 Gyr increments, since the initial 
models in the lower left. This plot shows why cooling takes longer with the new grids (in black). 
For Uranus, for instance, at 4.25 R e , a much smaller T im indicates a longer time step dt between 
successive models, as seen from Eq. @J Uranus, with its larger T eq , has much slower evolution in 
the past few gigayears, such that small changes in the atmosphere grids can lead to large changes 
in the evolution ages. 
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Fig. 1 5. — Luminosity of the giant planets as a function of age. Observed values dPearl & Conrath 
199 If) are shown as black dots. The values for Uranus and Neptune nearly overlap. Estimates 
at young ages are highly uncertain, and depend strongly on the details of the formation process, 
especially for Uranus and Neptune. Note that Uranus and Neptune cross at age of 1 Gyr, and 
Neptune becomes cooler. 



Table 1 . Jupiter grid at 1 and 0.7 L, 



gravity 



13.5 



7k 1.0 


71o, 1.0 


7k, 0.7 


Tie, 0.7 


Tim 


1201.6 


3471.1 


1201.6 


3471.1 


1200 


904.7 


3124.9 


904.7 


3124.9 


900 


751.3 


2894.2 


751.2 


2894.2 


750 


600.6 


2582.3 


600.5 


2582.4 


600 










450 










316 










224 










168 










133 










112 










100 










89 


1201.2 


3359.3 


1202.6 


3360.3 


1200 


902.5 


2984.3 


902.6 


2984.3 


900 


750.8 


2730.3 


750.7 


2730.2 


750 


600.4 


2391.7 


600.4 


2391.6 


600 


450.6 


1966.2 


450.4 


1966.1 


450 










316 










224 










1 /Z O 

168 










133 










112 










100 










89 


1201.8 


3255.1 


1202.2 


3255.8 


1200 


901.5 


2850.3 


901.5 


2850.2 


900 


750.5 


2571.4 


750.5 


2571.4 


750 


600.4 


2222.1 


600.3 


2222.1 


600 


450.6 


1816.0 


450.4 


1815.9 


450 


317.4 


1312.2 


317.0 


1310.2 


316 


227.3 


842.6 


226.4 


838.0 


224 










168 










133 










112 










100 










89 


1201.5 


3165.6 


1201.8 


3166.0 


1200 
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Table 1 — Continued 



gravity 7^,1.0 T w , 1.0 T^, 0.7 T w , 0.7 7- 



901.5 


2731.4 


901.0 


2731.4 


900 


750.4 


2435.8 


750.4 


2435.8 


750 


600.3 


2086.1 


600.3 


2086.0 


600 


450.6 


1701.6 


450.4 


1701.6 


450 


317.4 


1233.5 


317.0 


1231.7 


316 


227.4 


791.8 


226.4 


787.4 


224 


175.3 


558.7 


173.2 


552.4 


168 


146.4 


451.4 


142.7 


439.9 


133 


- 


- 


- 


- 


112 


- 


- 


- 


- 


100 


- 
- 


- 
- 


- 

- 


- 
- 


89 
1200 


- 


- 


- 


- 


900 


750.3 


2330.5 


750.3 


2330.5 


750 


600.3 


1987.6 


600.2 


1987.6 


600 


450.6 


1617.7 


450.4 


1617.6 


450 


317.4 


1174.5 


317.0 


1172.6 


316 


227.4 


754.4 


226.4 


750.4 


224 


175.4 


533.5 


173.2 


527.5 


168 


146.4 


430.9 


142.8 


419.6 


133 


131.9 


375.8 


126.8 


360.0 


112 


125.0 


348.9 


118.9 


330.1 


100 


119.7 


328.0 


112.9 


306.9 


89 










1200 










900 










750 


600.3 


1918.6 


600.2 


1918.6 


600 


450.6 


1561.0 


450.4 


1560.9 


450 


317.4 


1135.3 


317.0 


1133.5 


316 


227.4 


729.5 


226.4 


725.3 


224 


175.4 


515.9 


173.3 


510.0 


168 


146.5 


417.1 


142.8 


406.5 


133 


131.9 


364.3 


126.8 


348.5 


112 


125.0 


337.9 


119.0 


320.0 


100 


119.8 


317.9 


113.0 


296.9 


89 



Table 1 — Continued 



gravity 7^,1.0 T w , 1.0 7k, 0.7 T w , 0.7 7- 



25.1 











1200 




_ 






900 










750 










600 




_ 






450 


- 


- 


- 


- 


316 


227.4 


716.0 


226.4 


712.0 


224 


175.4 


506.9 


173.3 


501.1 


168 


146.5 


409.7 


142.8 


399.2 


133 


132.0 


357.4 


126.9 


342.3 


112 


125.1 


331.9 


119.1 


314.3 


100 


119.9 


312.5 


113.1 


292.3 


89 










1200 










900 










750 










600 










450 










316 


227.7 


703.7 


226.6 


699.5 


224 


176.0 


499.4 


173.7 


493.1 


168 


146.5 


402.3 


142.8 


391.6 


133 


132.0 


350.9 


126.9 


336.1 


112 


125.1 


325.8 


119.1 


308.7 


100 


120.0 


307.2 


113.2 


286.8 


89 



Note. — Surface gravities are in m s~ 2 . A metallicity of 
[M/H]=+0.5 (~ 3 x solar) is assumed. "0.7" and "1.0" mean 
0.7 and 1.0 times the current Jovian incident flux. 



Table 2. Saturn grids at 1 and 0.7 L Q 
gravity 7^,1.0 r 10 , 1.0 T eB , 0.7 T 10 , 0.7 7- 



1.3 


652.8 


3196.1 


652.8 


3196.1 


650 




551.5 


3020.1 


551.4 


3020.1 


550 




450.7 


2750.7 


450.7 


2750.6 


450 




350.5 


2232.2 


350.4 


2231.6 


350 



275 

225 

175 

131 

104 

87 

78 

69 



651.3 


2993.9 


651.3 


2993.8 


650 


550.7 


2783.9 


550.7 


2783.9 


550 


450.4 


2485.2 


450.4 


2485.2 


450 


350.3 


2000.7 


350.3 


2000.1 


350 


275.5 


1500.1 


275.6 


1498.8 


275 


226.1 


1120.5 


225.8 


1118.4 


225 



175 

131 

104 

87 

78 

69 



650.7 


2831.4 


650.7 


2831.4 


650 


550.4 


2597.4 


550.4 


2597.4 


550 


450.3 


2285.2 


450.2 


2285.1 


450 


350.4 


1844.9 


350.3 


1844.4 


350 


275.6 


1376.2 


275.5 


1375.0 


275 


226.1 


1021.8 


225.8 


1019.7 


225 


177.2 


692.2 


176.5 


689.2 


175 


135.7 


490.5 


134.3 


484.9 


131 


112.2 


380.9 


109.9 


373.3 


104 
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Table 2 — Continued 



gravity T eB , 1.0 T w , 1.0 7^,0.7 7i , 0.7 7- 



10.0 











550 










450 


350.4 


1769.9 


350.3 


1769.4 


350 


275.6 


1318.1 


275.5 


1317.0 


275 


226.1 


976.6 


225.7 


974.6 


225 


177.1 


661.7 


176.6 


658.7 


175 


135.7 


478.8 


134.3 


463.4 


131 


112.2 


364.8 


110.0 


356.5 


104 


99.6 


311.5 


96.0 


302.4 


87 










78 










69 










650 










550 










450 










350 










275 


226.1 


944.6 


225.8 


942.4 


225 


177.2 


640.2 


176.6 


637.9 


175 


135.8 


453.6 


134.3 


448.4 


131 


112.3 


353.3 


110.0 


345.1 


104 


99.5 


301.0 


96.4 


291.7 


87 


93.9 


279.2 


89.9 


268.2 


78 


0".z. 




84. 1 




(SO 










650 










550 










450 










350 










275 










225 


177.2 


628.6 


176.6 


625.7 


175 


135.8 


445.4 


134.4 


440.0 


131 


112.3 


347.5 


110.0 


339.2 


104 


99.5 


296.9 


96.1 


286.5 


87 


94.0 


275.9 


90.3 


264.2 


78 


89.2 


257.0 


84.1 


242.3 


69 










650 










550 
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Table 2 — Continued 



gravity 7k, 1.0 r 10 , 1.0 r eff , 0.7 T 10 , 0.7 T, 











450 










350 


- 


- 


- 


- 


275 










225 


- 177.2 


617.4 


- 176.6 


-614.9 


175 


135.7 


437.1 


134.4 


432.0 


131 


112.4 


341.7 


110.2 


333.3 


104 


99.5 


291.6 


96.1 


282.4 


87 


94.0 


270.3 


90.1 


259.9 


78 


89.2 


252.8 


84.2 


238.3 


69 










650 










550 










450 










350 










275 










225 










175 


135.8 


432.3 


134.4 


427.4 


131 


112.5 


338.1 


110.2 


329.7 


104 


99.6 


288.3 


96.2 


279.2 


87 


94.0 


267.3 


90.1 


257.1 


78 


89.2 


250.2 


84.2 


235.9 


69 



Note. — Surface gravities are in m s" 2 . A metallicity of 
[M/H]=+1.0 (lOx solar) is assumed. "0.7" and "1.0" mean 
0.7 and 1.0 times the current Saturnian incident flux. 



Table 3. Uranus and Neptune Grids at 4 Level of Solar Flux 



gravity T eS , 0.12N ri ,0.12N ri,0.12N r eff , 1.0N T w , 1.0N T u 1.0N T ef{ , 1.0U T w , 1.0U 7i, 1.0U 7k, 1.8 U 7i , 1.8U 7i, 1.8U 7- 



7.9 



32.66 


88.64 


49.11 


48.71 


129.32 


66.00 


59.87 


164.80 


80.55 


69.54 


196.07 


94.22 


27 


37.37 


105.39 


56.31 


50.38 


139.77 


70.22 


60.80 


171.55 


83.42 


70.14 


200.74 


96.34 


34 


44.89 


134.43 


68.06 


54.11 


160.32 


78.67 


63.04 


186.10 


89.75 


71.63 


211.46 


101.29 


43 


55.11 


177.20 


85.86 


60.87 


193.49 


93.05 


67.67 


211.78 


101.44 


74.30 


229.80 


109.99 


54 


68.46 


231.89 


111.00 


71.96 


243.06 


116.47 


76.39 


253.96 


121.88 


81.78 


267.16 


128.53 


68 


86.23 


301.80 


146.36 


88.21 


309.59 


150.42 


90.90 


315.21 


153.36 


94.32 


321.08 


156.43 


86 


109.12 


396.53 


196.11 


109.93 


399.17 


197.50 


111.22 


404.77 


200.44 


113.40 


409.02 


202.66 


109 


137.31 


531.03 


265.95 


137.48 


532.84 


266.87 


138.41 


535.66 


268.32 


139.49 


539.71 


270.40 


137 


217.06 


1086.00 


551.01 


217.15 


1086.72 


551.39 


217.22 


1087.84 


551.98 


217.63 


1090.67 


553.46 


217 


32.63 


85.80 


47.82 


48.59 


125.33 


64.40 


59.69 


159.82 


78.46 


69.34 


190.39 


91.66 


27 


37.35 


101.83 


54.82 


50.28 


135.23 


68.38 


60.63 


166.29 


81.18 


69.94 


194.96 


93.71 


34 


44.87 


129.67 


66.14 


54.03 


154.67 


76.32 


62.88 


179.89 


87.03 


71.44 


204.83 


98.21 


1 43 


55.09 


170.88 


83.13 


60.81 


186.84 


90.08 


67.53 


204.74 


98.17 


74.25 


222.41 


106.46 


654 


68.46 


224.03 


107.23 


71.88 


234.56 


112.30 


76.29 


245.17 


117.51 


81.63 


258.32 


124.06 


1 68 


86.23 


291.63 


141.08 


88.18 


298.61 


144.70 


91.04 


302.92 


147.98 


94.31 


310.17 


150.72 


86 


109.12 


383.20 


189.11 


110.26 


386.88 


191.04 


111.22 


391.57 


193.51 


113.02 


396.00 


195.83 


109 


137.31 


512.66 


256.50 


137.66 


514.23 


257.31 


138.17 


517.05 


258.75 


139.18 


520.87 


260.72 


137 


217.05 


1048.68 


531.44 


217.13 


1048.72 


531.48 


217.27 


1049.94 


532.10 


217.64 


1052.41 


533.39 


217 


32.60 


83.23 


46.61 


48.47 


121.55 


62.88 


59.53 


154.87 


76.40 


69.15 


184.65 


89.11 


27 


37.32 


98.37 


53.36 


50.21 


130.57 


66.51 


60.48 


160.84 


78.89 


69.76 


188.84 


90.97 


34 


44.85 


124.96 


64.25 


53.95 


149.20 


74.06 


62.74 


173.88 


84.43 


71.26 


198.27 


95.22 


43 


55.08 


164.66 


80.49 


60.75 


180.04 


87.10 


67.41 


197.53 


94.88 


74.20 


215.52 


103.18 


54 


68.70 


217.54 


104.15 


71.83 


226.18 


108.26 


76.19 


236.65 


113.32 


81.54 


249.39 


119.60 


68 


86.39 


282.31 


136.28 


88.14 


287.47 


138.93 


90.73 


293.36 


141.97 


94.31 


301.78 


146.35 


86 


109.12 


369.37 


181.83 


110.17 


373.36 


183.93 


111.20 


377.69 


186.21 


113.26 


383.08 


189.04 


109 


137.06 


494.06 


246.90 


137.48 


495.74 


247.77 


138.16 


498.64 


249.27 


139.34 


501.67 


250.83 


137 


217.02 


1010.68 


511.62 


217.09 


1010.97 


511.78 


217.22 


1012.28 


512.46 


217.60 


1014.69 


513.71 


217 


32.58 


80.44 


45.29 


48.38 


117.50 


61.25 


59.39 


150.10 


74.43 


68.98 


179.06 


86.67 


27 



Table 3 — Continued 



gravity 



i e ff, U.1ZJN 


1 10, U. 121N 


1 1, U. IziN 


T 1 fVM 


T 1 f\M 
ilO, l.UiN 


1 l , 1 .ULN 


T 1 f\I T 
Jeff, l.UU 


r r i rvi t 

i\o, i.uu 


T 1 f\I T 
11, l.UU 


l e fi, l.oU 


T 1 CI T 

iio, l.oU 


T 1 Oil 

il, l.oU 




37.31 


95.08 


51.95 


50.16 


126.26 


64.78 


60.34 


155.69 


76.74 


69.59 


182.81 


88.30 


34 


44.83 


120.47 


62.44 


53.89 


143.95 


71.92 


62.61 


168.05 


81.93 


71.11 


192.05 


92.40 


43 


55.06 


158.62 


77.96 


60.71 


173.66 


84.33 


67.30 


190.77 


91.83 


74.43 


209.41 


100.32 


54 


68.68 


209.94 


100.58 


71.77 


218.12 


104.41 


76.11 


228.40 


109.32 


81.44 


240.94 


115.42 


68 


86.38 


272.52 


131.26 


87.81 


277.19 


133.65 


90.87 


282.23 


137.19 


94.62 


290.64 


140.56 


86 


109.11 


356.83 


175.23 


109.92 


359.09 


176.42 


111.21 


362.70 


178.23 


113.01 


368.99 


181.63 


109 


137.09 


476.76 


237.96 


137.59 


478.66 


238.94 


138.34 


481.08 


240.19 


139.35 


484.79 


242.11 


137 


217.02 


974.45 


492.82 


217.09 


975.07 


493.14 


217.27 


975.94 


493.59 


217.59 


978.74 


495.04 


217 


32.54 


77.98 


44.08 


48.32 


113.63 


59.68 


59.28 


145.16 


72.41 


68.84 


173.51 


84.27 


27 


37.28 


91.67 


50.46 


50.05 


121.79 


62.98 


60.22 


150.46 


74.58 


69.44 


176.97 


85.76 


34 


44.82 


116.07 


60.67 


53.85 


138.64 


69.76 


62.50 


162.17 


79.44 


70.96 


185.53 


89.50 


i 43 


55.05 


152.70 


75.50 


60.60 


167.19 


81.56 


67.20 


183.96 


88.81 


74.34 


202.37 


97.08 


-^54 


68.66 


202.10 


96.96 


71.72 


210.16 


100.68 


76.06 


220.01 


105.32 


81.38 


231.67 


110.90 


1 68 


86.37 


262.93 


126.39 


88.03 


267.95 


128.93 


90.55 


272.98 


131.49 


94.17 


279.81 


134.99 


86 


109.12 


343.60 


168.27 


110.16 


346.52 


169.80 


111.18 


351.48 


172.41 


113.23 


355.11 


174.32 


109 


137.10 


458.98 


228.73 


137.62 


460.71 


229.63 


138.16 


463.29 


230.97 


139.31 


466.40 


232.58 


137 


217.02 


938.07 


474.00 


217.09 


938.65 


474.30 


217.22 


939.88 


474.93 


217.59 


942.07 


476.06 


217 


32.53 


75.61 


42.88 


48.26 


109.92 


58.17 


59.17 


140.49 


70.52 


68.72 


168.23 


82.01 


27 


37.27 


88.63 


49.11 


49.98 


117.66 


61.31 


60.13 


145.40 


72.51 


69.31 


171.54 


83.42 


34 


44.83 


111.89 


58.97 


53.78 


133.81 


67.81 


62.41 


156.53 


77.09 


70.84 


179.42 


86.83 


43 


55.03 


146.97 


73.15 


60.60 


161.00 


78.95 


67.20 


176.78 


85.68 


74.22 


195.32 


93.88 


54 


68.45 


193.69 


93.14 


71.70 


202.18 


97.00 


76.04 


212.20 


101.63 


81.23 


224.10 


107.26 


68 


86.23 


252.90 


121.35 


87.95 


257.84 


123.82 


90.52 


263.30 


126.57 


94.33 


271.78 


130.88 


86 


109.12 


331.47 


161.88 


110.09 


334.51 


163.48 


111.28 


337.63 


165.13 


113.00 


342.70 


167.79 


109 


137.06 


442.53 


220.17 


137.48 


444.04 


220.96 


138.30 


446.57 


222.28 


139.29 


449.89 


224.01 


137 


217.02 


902.64 


455.74 


217.09 


903.60 


456.24 


217.26 


904.67 


456.79 


217.59 


906.90 


457.93 


217 



Note. — Surface gravities are in m s 2 . A metallicity of [M/H]=+1.5 (~ 30 x solar) is assumed. "0.12N" and "1.0N" mean 0.12 and 1.0 times the current Nept 
incident flux, while "1.0U" and "1.8U" mean 1.0 and 1.8 times the current Uranian incident flux. The temperatures at 1 and 10 bars are provided. 



I 
I 



